Determination of the lowest energy structure of Ag 8 from first-principles calculations 
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The ground-state electronic and structural properties, and the electronic excitations of the lowest 
energy isomers of the Agg cluster are calculated using density functional theory (DFT) and time- 
dependent DFT (TDDFT) in real time and real space scheme, respectively. The optical spectra 
provided by TDDFT predict that the D2d dodecahedron isomer is the structural minimum of Ags 
cluster. Indeed, it is borne out by the experimental findings. 



PACS numbers: 36.40.-c, 36.40.Mr 

In latter years, a new field has emerged from the 
understanding, control and manipulation of objects at 
nanoscale level (nano-objects). It is commonly known 
as nanoscience. This field involves physics, chemistry 
and even engineering, and addresses a huge number of 
important issues starting from basic science and end- 
ing in a large variety of technological applications 0. 
Among the nano-objects, the small clusters or nanoclus- 
ters play a very important role, since they are the bricks 
of nanoscience. Therefore, the study of small clusters 
deserves a special attention. In this respect, the steps 
to follow for a complete description of a cluster can be 
summarized in the following three questions: what is the 
lowest energy structure?, what is the effect of increasing 
or decreasing the temperature on the structural proper- 
ties of a cluster? and the last step deals with the kinetic 
effects in the formation of the nanocluster. Here, we are 
only concerned about the first question for the silver oc- 
tamer, leaving the other two questions open for future 
investigations. 

From the theoretical point of view, first-principles 
methods give an enormous advantage for understanding, 
projecting and inventing new materials that is reflected in 
the huge number of articles published in the field of mate- 
rials science. Likewise, density functional theory (DFT) 
has emerged as a new and promising tool for ab initio 
electronic structure calculations and gives valuable infor- 
mation about the geometry of nanoscale systems [2j but 
unfortunately it not always predict the correct structure 
of the cluster under consideration. In this regard, the 
silver octamer belongs to the group of the controversial 
systems for which the lowest energy structure is unre- 
solved by DFT. 

Recently, P. Radcliffe et al. [3| have proposed Agg 
clusters embedded in helium droplets as a suitable sys- 
tem for light amplification based on an optically accessi- 
ble long- living excited state (E*) and thereby, from the 
theoretical point of view, the determination of the struc- 
ture becomes a key point as the first step for identifying 
and controlling the levels that populate E*. Up to now, 
it has not been possible to make a reliable theoretical 
prediction of the most stable structure of Agg . A review 



of the literature reveals that there are two competing 
geometries in eight-atom clusters of s-electron elements, 
having and T)2d symmetry. In fact, different levels 
of theories favor different geometries: DFT in its local 
density approximation (LDA) multireference config- 
uration interaction method |5J] , a tight-binding approach 
, and the many-body perturbation theory-based calcu- 
lations give the T>2d geometrical shape as the lowest 
energy structure of Agg whereas the equation-of-motion 
coupled cluster method time-dependent DFT only at 
LDA level and molecular-dynamics simulations |Hj 
predict a T^ structure as the structural minimum. It was 
reported in Refs. 0, that the D 2 d geometry is favored 
energetically over T^ symmetry when explicit correlation 
treatments for 5s electrons are included, but since the 
calculated energy difference between T^ and T)2d isomers 
is very small, the predicted theoretical ordering is un- 
certain. One way of solving this vexing problem comes 
from the hand of the time-dependent density functional 
theory (TDDFT) that is a generalization of tradi- 
tional ground stationary state DFT to treat the dynamic 
response of the charge density to a time-dependent per- 
turbation. TDDFT is a powerful methodology towards 
the calculation of the optical spectra, and thereby gives 
access to excited-state information. 

In this brief report, we have calculated the optical re- 
sponse of the Td and T>2d geometries and they were com- 
pared both with each other and the experimental find- 
ings. The atomic positions were fully optimized with 
an all-electron DFT implementation at the generalized 
gradient approximation (GGA) level, representing an im- 
provement over other TDDFT studies in small silver clus- 
ters . This work relies on the combination of the tra- 
ditional DFT and its generalization to excited states, as 
a promising tool for elucidating structures that DFT by 
its own is unable to predict. Here, we demonstrate that 
the Y>2d structure is the structural minimum of Agg and 
the calculated spectra allow us to estimate the interac- 
tion of Agg with the surrounding helium or argon matrix 
presented in experimental observations. 

With the aim of elucidating the lowest energy struc- 
ture of Agg cluster, we have performed density functional 
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theory-based calculations consisting of a linear com- 
bination of Gaussian-type orbitals-Kohn-Sham-density- 
functional methodology (LCGTO-KSDFM) to obtain the 
structural and electronic ground-state properties [l2^ . 
and a TDDFT implementation to compute the electronic 
excitations For the former, all-electron calculations 
were carried out with DEMON-KS3P5 [H at GGA level 
to take the exchange-correlation (XC) effects into ac- 
count 14]. An orbital basis set of contraction pattern 
(633321/53211*/531+) was used in conjunction with the 
corresponding (5,5;5,5) auxil iary basis set for describ- 
ing the s-, p- and d-orbitals [15j- The grid for numer- 
ical evaluation of the XC terms had 128 radial shells 
of points and each shell had 26 angular points. Spuri- 
ous one-center contributions to the XC forces, typically 
found in systems with metal-metal bonds when using a 
nonlocal functional, are eliminated in a similar way as 
has been done in Ref. ^(|. Trial geometries were fully 
optimized without symmetry and geometry constraints 
for different multiplicities using the Broyden-Fletcher- 
Goldfarb-Shanno algorithm [J^- The multiplicities were 
ranged from 1 to 11 and in all reported structures the 
singlet state was favored energetically. During the opti- 
mization, the convergence criterion for the norm of the 
energy gradient was fixed to 10~ 4 a.u. while it was 10~ 7 
a.u. for the energy and 10 -6 a.u. for the charge density. 
For the latter, after inserting the atomic coordinates of 
the converged structures provided by DEMON-KS3P5 
all the dynamical quantities are computed by evolving 
the electronic wave functions in real time and real space 
[l3j |. The electron- ion interaction is described through 
the Hartwigsen-Goedecker-Hutter relativistic separable 
dual-space gaussian pseudopotentials and the XC 
effect were treated in the GGA, implemented via the 
Perdew-Burke-Ernzerhof functional |19| . The grid in real 
space to solve the Kohn-Sham equations consists in a sum 
of spheres around each atom of radius 5.5 A and a mesh 
spacing of 0.23 A. The time step for the propagation of 
the electronic orbitals was fixed to 0.0013 fs, which en- 
sures the stability of time-dependent propagation. An 
artificial electronic temperature of 10 K was included ac- 
cording to the Fermi-Dirac function used to distribute 
the electrons among the accessible states. 

After a review of the literature on silver clusters 
01111 110, we have decided to optimize, as a good 
candidate to the structural minimum of the octamer, the 
following isomers of Ag 8 : a D 2( j dodecahedron (D 2( i-DD), 
which can also be viewed as a distorted bicapped octa- 
hedron, a Td tetracapped tetrahedron (T^-TT) and a 
C s 1-pentagonal bipyramid (C S -PBP) in Fournier's no- 
tation 0. The main results (density of states (DOS), 
optimized structures, polarizabilitics, ground state en- 
ergies, . . . ) of the electronic structure calculations are 
collected in Table [Q and Fig. □ The LCGTO-KSDFM 
calculations clearly show that the C S -PBP geometry is 
energetically far from the lowest energy structure by an 



TABLE I: Ground-state energies relative to the most stable 
isomer (D2d-DD) and electronic structure properties of the 
DFT-optimized Ag8 cluster isomers. The Fermi level is de- 
noted by E f and A£ stands for the HOMO-LUMO gap. The 
mean static polarizability a and the polarizability anisotropy 
Aa were calculated under the influence of an external electric 
field of strength 0.0005 a.u.. 



Symmetry 


AEdft 


E/ 


AC 


ex. 


Aa 




(meV) 


(eV) 


(eV) 


(A 3 /atom) 


(A 3 /atom) 


E>2d 


0.00 


-4.302 


1.719 


6.32 


1.22 


Td 


6.14 


-4.574 


2.335 


6.46 


0.01 


Cs 


181.25 


-4.048 


1.327 


6.45 


1.71 



amount of 181.25 meV. Therefore, we will concentrate 
our attention in Y>id and T^ structures. 

Despite the fact that the lowest energy isomer corre- 
sponds to a D2d symmetry, it should be noted that the 
ground-state energy difference between D2<j-DD and T^- 
TT isomers is very small (AE]j 2d ^T d =6.12 meV) com- 
pared to 0.19 eV which is the averaged energy difference 
between the ground-state structure and the second sta- 
ble structure of Ag„ (2<n<12) 0. Furthermore, the po- 
larizabilities and the HOMO-LUMO gaps (HLg) do not 
offer a clear picture for elucidating the structural mini- 
mum of the octamer. That is, it is well known that the 
Ags cluster is a closed-shell system and it was demon- 
strated experimentally and theoretically that the closure 
of electronic shell manifest itself in particularly large HLg 
( see jl| and references therein), consequently the HLg 
reported in Table [i] arc on the side of stabilization of 
the T^-TT isomer. As far as the reported polarizabil- 
ities are concerned, in molecular electronic distribution 
studies under the influence of an external electric field, 
the relevant quantities are the mean static polarizability 
a = (E?=i a ")/3 and the polarizability anisotropy Aa, 
defined as: 
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where otij = d(ji e )i/dEj is the ij-component of the po- 
larizability tensor under the action of an external elec- 
tric field Ej. It is not a common procedure to express 
the polarizability anisotropy such as it was defined in 
Eq. QJ. The commonly- used definition omits the second 
term (6 JT^ ' J=1 2 -i<j a %) an( ^ thus ne gl ec ts the important 
influence that the off-diagonal elements of the second- 
rank polarizability tensor play in the symmetry consid- 
erations of the electric charge distribution [2(| . The mean 
static polarizabilities, reported in Tabled are quite sim- 
ilar to each other showing in average that the electron 
charge is nearly equally distributed among the three iso- 
mers. However, only the polarizability anisotropy of the 
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FIG. 1: (color online). Energy levels, partial and total density of states and shapes of the delocalized molecular orbitals for 
the higher-lying occupied and lower-lying unoccupied levels of Agg isomers: a) Td-TT and b) D2d-DD. The orbitals below the 
Fermi level are double occupied; for these levels only one molecular orbital is presented. The solid vertical line represents the 
Fermi level. 



Trf-TT isomer is clearly reduced. This result tends to 
stabilize the T^ symmetry over D 2 d because the less po- 
larizability anisotropy is, the more spherically symmetric 
charge distribution is and the latter condition is favored 
by a closed-shell system like the silver octamer. Indeed, 
the delocalized molecular orbitals for the higher-lying oc- 
cupied levels of the D 2 d and T^ geometries, presented in 
Fig- IH exhibit a hybridization of the atomic 5s levels with 
the 5p levels leading to a nearly spherical shape whereas 
for the lower-lying unoccupied levels, the spherical sym- 
metry becomes less important. The proximity in energy 
and the corresponding superposition of the spatial distri- 
bution of the higher-lying occupied molecular orbitals of 
the Trf-TT isomer with respect to the ones of the D 2( rDD 
isomer, contribute to a spherical symmetrization of the 
Trf-TT charge distribution. 

The balance of the aforementioned contrary tenden- 
cies does not allow a reliable prediction of the structure, 
such as it was discussed in literature 0, EJ 0, H, Ell- 
The TDDFT calculations using the structural parame- 
ters provided by the LCGTO-KSDFM as starting point 
can shed some light on the better understanding of this 
vexing controversy on the structure of Ags. In this re- 
spect, as it is shown in Fig. the calculated spectrum 
for D 2( j symmetry is in excellent agreement with the res- 
onant two-photon ionization spectrum reported by F. 
Federmann et al. [2l| and the excitation spectrum re- 



ported by C. Felix et al. [13], whereas the T rf -TT cal- 
culated spectrum only shows one resonant peak and it is 
about 0.31 eV blue-shifted with respect to the experimen- 
tal measurements. Thus, the LCGTO-KSDFM predicted 
structure is confirmed by the TDDFT calculations of the 
optical response when it is compared to the experimental 
evidences. 

Some of the relevant states involved in the transitions 
that populate the peaks of Fig. |21 are provided through 
the calculated DOS depicted in Fig. ^ For the D 2 d iso- 
mer, the eigenvalues of the HOMO and HOMO-1 states 
are both close together in energy and the HOMO-2 state 
is 0.7 eV further down meanwhile for the T^-TT isomer, 
this three states were grouped together in a window en- 
ergy of only 3 meV. It is worthwhile to mention here that 
the Td-TT isomer has high symmetry, so many states will 
be degenerate. However, because of small numerical er- 
rors, it is quite possible that states that should strictly 
be degenerate will show as being within very small en- 
ergy windows. In the energy range displayed in Fig. |3 
for D 2 d symmetry, the TDDFT calculations predict four 
electronic transitions starting from the twofold degener- 
ate HOMO and HOMO-1 states to two excited states sep- 
arated 0.07 eV in energy and it gives rise to two peaks 
because of the nearness in energy of the HOMO-1 and 
HOMO states. In the case of the T^ symmetry, six tran- 
sitions are predicted giving rise to only one peak because 
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FIG. 2: Comparison between two experimental recorded spec- 
tra and the calculated spectra for both (solid line) r>2d- 
DD and (dashed line) T^-TT isomers at a temperature of 
10 K. Open circles correspond to the resonant two-photon- 
ionization (R2PI) spectroscopy on Ags clusters in He droplets 
PH while the solid triangles are for the excitation spectrum 
of Ag8 excited with monochromatic Xe light in an Ar matrix 
(22^. The D2d-DD isomer spectrum is in excellent agreement 
with the experiment whereas the TVTT isomer spectrum is 
around 0.31 eV blue-shifted compared to the R2PI spectrum. 

the two excited states that electronic transitions populate 
are only 1 meV separated, as commented above. 

Some consequences can be extracted going further in 
the analysis of the calculated spectrum for D2<j-DD iso- 
mer when it is compared with the experimental references 
reported in Fig. [3 On one hand, for the case of R2PI 
experiment we attribute the slight difference in peak po- 
sition (~ 3 meV) to the helium environment through the 
formation of electron bubble states that significantly blue 
shift the transition 23]. As already mentioned above, 
our TDDFT calculations also confirm the authors' feel- 
ing of Ref. [2l| that the asymmetry of the peak involves 
more than one transition. On the other hand, the com- 
parison between the excitation spectrum in Ar matrix 
|22| and the calculated T>2d isomer spectrum allow us to 
measure the interaction of the Ar matrix with the Ags 
cluster. Thus, the shift of energy probably due to Ar ma- 
trix effects is estimated to be about 25 meV that is great 
enough compared to the slight difference in peak position 
(~ 3 meV) attributed to the helium surrounding the Agg 
clusters in R2PI experiment. Consequently, the use of 
liquid helium droplets as a spectroscopic matrix has the 
advantage over argon matrix of providing an environment 
more suitable for study the electronic excitations of small 
and free silver clusters. 

In conclusion, we have shown that a combination of 
a LCGTO-KSDFM and TDDFT approach is able to re- 
produce the measured optical response of the silver oc- 
tamer and allow us to elucidate its lowest energy struc- 



ture below 10 K. Our calculation thus confirms that the 
structural minimum of Ags is the D2d-DD isomer, whose 
geometrical structure is depicted in Fig.^ The TDDFT 
calculations have provided a number of electronic transi- 
tions involving the resonant peaks showed in Fig. [21 and 
demonstrate that the R2PI experiment is a good tech- 
nique to measure experimentally the electronic excita- 
tions of bare silver clusters because it is less aggressive 
than, for example, the experiments that consider argon 
as the spectroscopic matrix. 
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